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ABSTRACT 

Observations of X-ray cavities formed by powerful jets from AGN in galaxy clus- 
ter cores are widely used to estimate the energy output of the AGN. Using methods 
commonly applied to observations of clusters, we conduct synthetic X-ray observations 
of 3D MHD simulated jet-ICM interactions to test the reliability of measuring X-ray 
cavity power. These measurements are derived from empirical estimates of the enthalpy 
content of the cavities and their implicit ages. We explore how such physical factors 
as jet intermittency and observational conditions such as orientation of the jets with 
respect to the line of sight impact the reliability of observational measurements of cavity 
enthalpy and age. An estimate of the errors in these quantities can be made by directly 
comparing "observationally" derived values with "actual" values from the simulations. 
In our tests, cavity enthalpy derived from observations was typically within a factor of 
two of the simulation values. Cavity age and, therefore, cavity power are sensitive to the 
accuracy of the estimated inclination angle of the jets. Cavity age and power estimates 
within a factor of two of the actual values are possible given an accurate inclination 
angle. 

Subject headings: galaxies: jets - galaxies: clusters: general - methods: numerical - 
X-rays: galaxies: clusters - magnetohydrodynamics (MHD) 



Introduction 



X-ray images of giant cavities in galaxy clusters associated with powerful jets from central 
active galactic nuclei (AGN) suggest that AGN may play an important role in the energetics of 
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The estimated minimum energy required to produce cavities is often in the range of 10 55 to 10 ( 



60 



erg (Birzan et al. 2004). Observations of ICMs have shown the existence of a temperature floor of 



approximately 2 keV {e.g., Peterson et a/.|[2002 ). The lack of gas below this temperature, contrary 
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to expectations of a classical "cooling flow" (Fabian 1994), is historically known as the "cooling 



problem" . Evidently, the energy required to suppress this cooling of the ICM below 2 keV is on the 



same order as energy in the X-ray cavities (McNamara & Nulsen 2007). As a result, one popular 



hypothesis that has emerged to solve the cooling problem is that energy injected into the ICM 
by AGN will quench cooling and subsequent star formation in the central cluster galaxy. Several 



numerical studies such as Briiggen et al. (2005) and Sijacki & Springel (2006) have supported this 
hypothesis. 

X-ray cavity systems are evidently formed when low density, hot plasma originating from the 
AGN inflates a bubble in the ICM. The low density plasma produces a decrement in the line 



of sight intensity through the cavity from the normal ICM X-ray emission (Clarke et al. 1997). 
The cavities produce roughly elliptical brightness depressions ~20% to 40% below the surrounding 
regions ( McNamara &: Nulsen|[2~007 ) . A few dozen such cavity systems are known (see Dong et al 



2010 Rafferty et a/.|2006 . and references therein). Some cavities are filled with radio emission from 
relativistic particles and are typically found in pairs with an AGN in the cluster center between 
them. Other cavities devoid of radio emission above 1.4 GHz are referred to as radio ghosts; (e.g., 
Birzan et a l. 2008). The presence of multiple cavity pairs in some cases suggests a series of outbursts 
from the AGN. Hydra A, for example, contains several attached cavities filling at least 10% of the 



cluster volume within 300 kpc of the cluster center (Wise et al. 2007). Work done by Morsony et 



al. (2010), however, suggests that the presence of multiple cavities may be the result of the motion 



of a dynamic ICM. The size of cavities varies greatly from 1 kpc in diameter for M87 (Young et al. 



2002), for example, to over 200 kpc in diameter in Hydra A (Wise et al. 2007). 



X-ray cavities are likely to be long lived structures, remaining intact for over 100 Myr in Hydra 



A, for example (Nulsen et al. 2005). On the other hand, simple hydrodynamic analyses suggest 



cavities filled with light gas should be unstable to Rayleigh Taylor (RT) and Kelvin Helmholtz (KH) 
instabilities as they form and rise in the cluster. Several numerical studies have been performed, 



which include additional physics to stabilize the bubbles. Jones & DeYoung (2005), for example 



carried out 2D calculations of bubbles with magnetic fields finding that the fields suppress instabil- 
ities. Reynolds et al. (2005) demonstrated the stabilizing effect of a Braginskii viscosity mitigated 



by Coulomb collisions. Briiggen et al. (2009) included a model for RT driven sub-grid turbulence 



in 3D hydrodynamic simulations of bubbles. Their results show that turbulence can also prevent 
the break up of bubbles as a by-product of resulting ICM entrainment. 

Surveys of cavity systems and their energy of formation requirements have found that nearly 
half of the studied cavities show evidence for sufficient power to suppress cooling for short periods 



of time in their host clusters (Birzan et al. 2004 Rafferty et al. 2006) if the energy in the cavities 



becomes distributed in the ICM. There are, however, significant uncertainties in the determination 
of the cavity energy contents and associated time scales. The energy content, generally assumed 
to be measurable in terms of the supporting pressure of the cavity, requires, for instance, accu- 
rate measurement of the cavity pressure and also its volume. The cavity volume, V, is generally 
estimated from circles or ellipses fit by eye to the cluster X-ray surface brightness distribution. 
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That two dimensional projection is then converted into a three dimensional ellipsoid by revolution 
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to inflate a cavity containing internal energy, E, is taken to be the enthalpy content of the cavity, 
E + PV, based on the assumption that the cavity expanded subsonically in the ICM at constant 
pressure, P. The timescale needed to innate the cavity and to measure the associated cavity power 
is usually estimated from buoyancy or characteristic sound crossing times arguments. 

A substantial amount of effort has gone into numerical studies of outflows from AGN and 
the creation of bubbles containing hot AGN generated plasma. To make direct comparisons with 
observations, however, realistic synthetic observations from these calculations are needed. Several 
authors have applied this approach with various models of jets and bubbles in a cluster (e.g., 



Briiggen et al. 2005 Diehl et al. 2007 Briiggen et al. 2009 Morsony et al. 2010). The synthetic 



observations of magnetically dominated cavities by Diehl et al. ( 2007 1 were able to produce several 
of the observed characteristics seen in real observations including bright rims commonly found 
outlining cavities ( McNamara fe Nulsen||2007 ). Briiggen et al. (2009) were able to show consistency 
between their measurement of PV from synthetic observations of bubbles with sub-grid turbulence 



and a sample of observed cavities by Diehl et al. (2008). 



Studies such as Dong et al. (2010) and Enfilin & Heinz (2002) have tested the efficiency of 



detecting cavity systems from X-ray observation, but to our knowledge, a detailed assessment 
of the reliability of the observational techniques used to determine cavity enthalpy has not been 
performed. Synthetic observations of the complex interactions involved in the formation of X- 
ray cavities provide a powerful test of these methods. The primary goal of this paper is to test 
common observational techniques for determining cavity energetics. We employed a pair of 3D 
magnetohydrodynamic (MHD) simulations of jets in realistic cluster environments presented in 



O'Neill & Jones (2010) (henceforth, OJ10). These simulations were post-processed to yield synthetic 
X-ray observations. Section[2]describes the models and numerical methods. Section[3]and ^4]present 
the observations and analysis, while §[5] lists the conclusions of this work. In the analysis we have 
used H = 72 km s _1 Mpc -1 , Sl M = 0.3, and ft A = 0.7. 



2. Simulation Details 



The simulations presented, described by OJ10, were computed on a 3D Cartesian grid using 



a 2nd order total variation diminishing (TVD) non-relativistic MHD code described by Ryu & 
Jones (1995) and Ryu et al. (1998). A gamma-law equation of state was assumed with 7 = 5/3; 



radiative cooling was negligible for the conditions of these simulations and was therefore ignored. 
Computational details are provided in OJ10. We provide here only an outline as needed to evaluate 
the present work. The physical extent of the computational grid was x = 600 kpc, y = z = 480 
kpc. Each computational zone represented one cubic kiloparsec with Ax = Ay = Az = 1 kpc. Two 
oppositely directed jets were centered within the grid and aligned with the x-axis. A passive tracer, 
Cj e t, was advected with the flow to identify jet material from ambient material. Two different jet 



models are utilized in the present discussion; 1) a so-called relic (RE) model in which quasi-steady 
jets were on for 26 Myr then turned off and, 2) an intermittent (113) model in which the jet power 
cycled on and off at 13 Myr intervals throughout the simulation. 

2.1. Bi-directed Jet Properties 

Both simulations featured bi-directed jets that had an internal Mach 3 speed at full power, 
corresponding to a physical speed Vj e t = 0.10c. These jets originated from a cylindrical region r Je t = 
3 kpc in radius and lj e t = 12 kpc in length centered in the grid. The gas injected at the jet origin was 
less dense than the ambient gas by a factor of one hundred and was initially in pressure equilibrium 
with its local surroundings. Temporal variation in the jet was controlled by an exponential ramp 
in density, pressure and momentum density over 1.64 Myr for 113 and 0.65 Myr for RE. Physical 
conditions inside the jet source region were relaxed to a volume average from a sphere surrounding 
the jet origin as the jet turned off, then evolved back from instantaneous volume averages for the 
local medium to the desired jet conditions as jet power resumed. The combined power from both 
jets at peak was L = 1.2 x 10 46 erg s . Small magnetic and gravitational energy contributions 
to the jet energy flux were ignored in defining the jet power. Those energy terms were, however, 
followed explicitly in the simulations and accounted for in energy exchanges between the jets and 
their surroundings. 

The magnetic field launched from the jet was purely toroidal, B§ = Bo(r/rj et ) inside a jet core 
region, with (3 = Pj e t/(B 2 /8ir) ~ 100, on the perimeter of the jet core. There was a thin 'sheath" 
surrounding the core, through which all the jet properties, including the magnetic field transitioned 
to local ICM conditions. 

2.2. Cluster Environment 

The simulation cluster environments in OJ10 were designed to mimic a realistic, relaxed cluster. 
Gravitational potential and density profiles were selected to yield a temperature profile typical of 
clusters in hydrostatic equilibrium. A tangled ambient magnetic field with a characteristic coherence 
length typical of observed clusters was chosen to break symmetry over the grid. The local ICM 
magnetic pressure averaged to about 1% of the gas pressure, although that ratio fluctuated by large 
factors over the volume. 



The NFW (Navarro et al. 



1997) dark matter density distribution 
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was used to generate the gravitational acceleration 



In 1 + 



rdn 



r + r dn 



(2) 



where r^ m = 400 kpc and p s ~ 4.3 x 10 g cm". This gave a virial mass of M v = 5 x 10 M G 







for a virial radius of 2 Mpc, which was within a factor of a few of the Perseus Cluster (e.g., Ettori 



et al. 1998) 



The gas density of the cluster was initialized with a density distribution 



Pa{r) = Po 
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given by OJ10, where /i = 0.9, f% = 0.1, r c i = 50 kpc, r C 2 = 200 kpc and /3 = 0.7. The density scale 
was po = 8.33 x 10 -26 g cm -3 . The pressure was determined by hydrostatic equilibrium, yielding a 
temperature profile resembling typical clusters (cf OJ10). The central pressure was Po = 4 x 10~ 10 
dyne cm -2 giving a sound speed Co = 895 cm s _1 in the cluster core for a 7 = 5/3 gas. 

On top of the initialized hydrostatic equilibrium in the ICM, a Kolmogorov spectrum of density 
fluctuations was imposed with a maximum local amplitude, ±0.10p a (r), as described by OJ10. 

The initially tangled and divergence-free cluster magnetic field was given by OJ10 as 

t = B e § + (4) 



where the components are 



Ba 



F\ (r) • m 



svno cosymc 
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with m = n = 3. F\(r) and i ? 2( r ) are functions designed to keep an approximately constant /3 
atmosphere with fluctuations that vary over scales of a few tens of kpc. The scale of the fields 
maintains a (3 ~ 100 on average over the cluster volume. The maximum magnitude of the field is 
~ 10//G. 



2.3. Relativistic "Cosmic Ray" Electrons 



The simulations included a population of relativistic Cosmic ray electrons (CRs) passively 
advected with the MHD quantities. The numerical details of the CR transport are given in Jones 
7taL\ fll999"l) , |Tregillis et at] (|200T|) and |Tregillis, Jones, Ryu| ([2004]) . A small, fixed fraction of the 
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thermal electron flux through shocks was injected into the CR electrons population and subjected 
to first order Fermi acceleration according to the standard test-particle theory. Downstream of 
shocks the CRs were also subject to adiabatic and synchrotron/inverse Compton radiative energy 
changes. The nominal CR pressure, which was neglected, was generally less than 1% of the gas 
pressure. CRs with Lorentz factors from 7 = 10 to 7 ~ 1.6 x 10 5 were tracked as a piecewise power 
law distribution. 

The inclusion of CRs allowed us to calculate inverse Compton and synchrotron emissions in 
a self-consistent manner. A separate analysis paper will include detailed consideration of radio 
synchrotron emission and high energy non-thermal X-ray (> 10 keV) emission. Only X-rays below 
10 keV are discussed in this paper. Those are entirely dominated in our computations by thermal 
emissions from the cluster ICM, although inverse Compton emissions are included. 



3. Synthetic X-ray Observations 

The physical quantities evolved through MHD simulations of radio jets provide unparalleled 
intuition into the complex dynamics of MHD flows, but it has not always been intuitive to relate 
these quantities to observation. To make this connection and address questions raised from obser- 
vations, emission processes in these simulations must be properly calculated and converted into a 
synthetic observation. 

The approach used here to model synthetic X-ray observations was based on |Tregillis et al.\ 



(2002). Observations were computed for an assumed cluster redshift, z = 0.0594, (Dl = 240 Mpc) 



corresponding approximately to the Hydra cluster ( Wise et al.|2007 ) . In order to understand better 



the influence of projection effects, we carried out the synthetic observations at three representative 

OO o 

angles, i = 80 ,45 , and 30 , between the jet axis and the line of sight. Two emission mechanisms 
were included; thermal bremsstrahlung and inverse Compton scattering off of CMB photons. In 
each zone of the computational grid we calculated emissivities based on local properties of the 
thermal or CR electron population, then corrected them for the cluster redshift. 

Thermal bremsstrahlung or free-free emissivity was computed as 

j Ulocal = 5.4 x lO- 39 5// (^ cai ,T e )Z^ e -^/^ ergcm^s^sr- 1 , (7) 

where uiocai = v obs (1 + z )- The free-free Gaunt factor, gtt, was computed by interpolation from 



the values calculated for plasma with typical ICM properties in Table 1 of Nozawa et al. (1998). 
We assume a fully ionized Z{ = 1 hydrogen gas with an ideal gas equation of state where the 
average temperature per zone was T e = Tj = T(keV) = /i-Pm#/(1.602 x 10 _9 p) with n = 1/2 
and P and p in cgs units. The numerical resolution of discontinuities in the simulations was a few 
zones. Consequently, the contact discontinuity between AGN (jet) and ICM plasmas was a few 
zones of moderately high density, very high temperature gas. These transition regions were artificial 
and should not in the absence of some equivalent, real viscous mixing, contribute to line of sight 
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intensities. To reduce this artifact, any zone with Cj e t > 0.01 (partially AGN plasma) had jv local for 
thermal bremsstrahlung set to zero. Equation [7] was integrated numerically over a given range of 
frequencies to simulate finite band widths of real instruments. This paper focuses on energy ranges 
accessible to observatories typified by Chandra. At those energies the inverse Compton emission in 
these simulations is negligible. Consequently, we omit details of their computation. 

Assessing projection effects is critical to comparisons of these synthetic observations with real 
observations. We developed a parallelized ray casting engine that allows the user to define an 
arbitrary orientation and resolution for the output images. A ray was cast normal to the image 
plane through the appropriately aligned grid of emissivities. Tri-linear interpolation was used at 
regular intervals along the ray and summed to give the total intensity along the line of sight, 
assuming an optically thin medium. Finally, intensities were converted into fluxes per pixel by 
multiplying the line-of-sight intensity by the solid angle of an image pixel. The image resolution 
was set to 1 arc sec, which matched the simulation 1 kpc physical resolution at the selected 240 
Mpc source distance. 



3.1. Relic (RE) Observations 

Synthetic X-ray observations of the RE, relic simulation are shown for several times and pro- 
jection orientations in Figure [TJ Following a common practice designed to highlight AGN-blown 
cavities, the computed brightness distribution of the ICM outside of the identified cavities was fit 
with a double /3-profile (§ |4.1| ) and then divided out to accentuate the X-ray cavities. The double 
/3-profile was determined independently for each synthetic X-ray observation. Figure [T] shows the 
synthetic X-ray observations in a 1.5-2.5 keV band divided by the best-fit double /3— profile in each 
instance. Time evolution of the system is displayed from left to right. At the earliest time shown, 
26.3 Myr, the jets had just turned off. Each row corresponds to a different orientation, with the 
inclination angle of the jet with respect to the line of sight decreasing from top to bottom. There 
are several notable features in each observation. Cavities are seen as brightness decrements from 
the surrounding emission. A pair of cavities is seen at 26.3 and 52.5 Myr at large inclination but 
appear to merge into a single cavity at small inclination. Presumably X-ray emission from a central 
galaxy would prevent the two cavities from appearing as a single cavity. Our simulations, however, 
do not include emission from gas bound specifically to the central galaxy. All inclinations show 
a pair of cavities at 157.5 Myr. The contrast in brightness of the cavity to the surrounding gas 
diminishes with both distance from cluster center and decreasing inclination. A detailed discussion 



on these trends can be found in Enfilin & Heinz (2002). The bow shock from the jets is seen in all 



observations. At early times the bow shock appears as a bright rim surrounding the cavities. At 
later times the bow shock has moved far from the edge of the cavities and no longer appears as a 
bright rim. 
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3.2. Intermittent (113) Observations 



Figure [2] shows a similar set of observations to those in Figure [TJ but for the 113, intermittent 
jet simulation. Several new features are seen with the introduction of jet intermittency. Late 
times at every inclination reveal "ripples" between the cavities and the bow shock. These features 
correspond to sound waves generated at the cavity walls during periods of jet activity. Similar 



"ripples" have been seen in observations of the Perseus cluster (Fabian et al. (2003)). A second 



related distinction from the RE observations is the appearance of bright rims outlining the cavities 
at every epoch. At smaller inclination angles the bright rims resemble the "arms" seen on smaller 



scales in NGC 4636 (Baldi et al. 2009) 



4. Cavity Measurements 

In the following analysis we attempted to apply common techniques for extracting two fun- 
damental parameters for each cavity detected in each observation throughout the elapsed time of 
each jet model. Every epoch for a given simulation represents a separate test for measuring both 
cavity enthalpy and cavity age. Following this time evolution allows us to detect biases and trends 
in the quality of the measurements. For the remainder of this paper we refer to values measured 
directly from the simulation data as the "actual" values, while values measured from the synthetic 
observations are referred to as "observed" values. We report the fractional error on a measured 
quantity x as e x = (x observed - x actual) /x actual for the remainder of this paper. 



4.1. Enthalpy 

The minimum energy required to produce a cavity is generally estimated as the total thermal 
energy in the cavity and the work done inflating the cavity slowly at constant pressure; that is, the 
enthalpy in the cavity, H = Utherm + PV ~ several x PV. In particular, if the adiabatic index of 
the cavity plasma is 7 C; 

H = -^^PV. (8) 

1c - 1 

For a gas with 7 C = 5/3, applicable to our simulations, this gives H = (5/2)PV. Estimation 
of cavity enthalpy, under the assumption that the cavity was inflated at its current location, re- 
quires knowledge of both the cavity volume and surrounding gas pressure. Since the AGN activity 
disturbed large volumes of the ICM it is not straightforward to determine either its pressure distri- 
bution or, for that matter, the volume occupied by the AGN generated cavity. A common strategy 
to resolve these two problems involves fitting a simple, symmetric brightness profile to regions of 
emission that seem not to include cavity structures. That profile can then be used to obtain es- 
timates for the average radial ICM properties. There are several variations of this strategy (e.g., 
Wise et aT]|2004 Birzan et aZ.||2004 ) . Our goal was not to determine the best strategy but to use a 
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common approach as an example. We followed a procedure similar to Wise et al. (2004) and Xue 



& Wu (2000) to extract pressure and Birzan et al. (2004) to extract volume from each observation 



For this exercise we used a double /3— profile profile {e.g., Ikebe et al. 1996) of the form 

l/2-3/3i 



Sx(r P ) = So Si 



01 



1 + 



02 



1 + 



R C 2 



-i l/2-3/3 2 N 



(9) 



where r„ is the projected distance from cluster center, to model the brightness distribution of the 



X-ray emitting ICM. Xue k, Wu (2000) discuss the benefits of using this profile as opposed to a 



single (3— profile. The profile was fit independently to each 1.5-2.5 keV synthetic X-ray observation 
of the RE and 113 simulations. The synthetic images were divided into annular bins, each ~ 1 arc 
sec in width. To remove any effects of the X-ray cavities in characterizing the brightness profile of 
the cluster plasma, a set of ellipses was chosen that best fit each cavity by eye. Any pixels within 
these ellipses were excluded from the annular bins. The average flux from the remaining pixels was 
used to define an azimuthally averaged brightness profile that was fit with the double /3-profile. 
Refer to Appendix [A] for details regarding the fitting procedure. Figure [3] shows example double 
j3— profile fits for observations of both models at an inclination of i = 45°. The best fit profiles 
resulted in 0.5 < Pi < 1.5, 0.9 < p 2 < 1.8 with typical values R C i ~ 50 kpc, R C 2 ~ 200 kpc, 
Soi ~ 0.8, and S02 ~ 0.2. The undisturbed cluster parameters were Pi = 0.7, Pi = 1, Rci = 55 
kpc, Rc2 = 260 kpc, Sqi = 0.9, and S02 = 0.1. 



4-1.1 . Cluster Temperature Profile 

The ICM temperature, Tjcm, at a given projected radius, r p , was determined from the ratio 
of fluxes in two bands; 1.5-2.5 keV and 9.5-10.5 keV. In particular, the equation 

Sx,i,-2,(r P ) = mh TiCM)e-»»l T >™<h, 



was solved for Tjcm using the aforementioned double /3-profile fits. Following Wise et al. (2004), we 
assumed that the two components of the double /3-profile corresponded to two phases of the ICM 
with temperatures Tjcm,\ for the inner component and Tjcm,2 for the outer component. Tjcm,i 



was taken to be the minimum and Tjcm,2 t ne maximum temperatures found using Equation 10 
The projected radius for the transition from Tjcm,i to Tjqm,2 was chosen to be the average of Rqi 
and Rc2- Figure [4] shows a comparison between the actual azimuthally averaged temperature profile 
as a function of physical radius from the RE initial conditions and the two component projected 
profile. Note that Tjcm,i mostly exceeded the actual inner core temperatures, since hotter gas 
along the line of sight contaminated Sx at small r p . 
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4-1-2. Cluster Electron Density Profile 



The radial thermal electron density profile of component i = 1,2 was obtained by inverting 
equation ^ , following the derivation of Xue & Wu ( 2000 ) ; 



n 2 ei (r p = 0) 



4tt 1 /2 



a(T ICM ,i)giVe \F(3f3i - 1/2) J \R Ci 
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i j • 



where 



a{Tj C M,: 



2 4 e b / 2tt 1.602 x 10- y Tf CM; 



3m e hc 2 

l/=(l+z)2.5keV/h 



3m P c 2 



1/2 



u=(l+z)1.5keV/h 



9ff 



^,T IC M,i)e- hu/TlCM dv, 



and 



3_ _ , RaSoj9i f TlCM 
hj RcjSoidj \TlCM,j 



A 1/2 



r(3/3,)r(3A - 1/2) 

r(3A)r(3/3, - 1/2) 
j = 1,2 and j ^ i. 



(11) 



(12) 
(13) 



(14) 



The values for Soi, Rci, and /3, were the best fit values for each component from the 1.5-2.5 keV 
observation. For simplicity, we assumed pure hydrogen. The electron weight, [i e = 2/(1 +X), where 
X, the hydrogen mass fraction, was therefore unity. The total electron density at a projected radius 
r p was determined by 



n e(r p ) = ^n ei (r p ) = n e (r p = 0) ^n«(r p = 0) 



i=l 



i=l 



1 + 



-3ft 



1/2 



(15) 



Figures [5] and [6] show example observed electron density profiles determined by Equation 15 
compared to the actual azimuthally averaged electron density from the RE simulation from obser- 

o 

vations at i = 45 . The data for the actual density were generated considering only computational 
zones with Cj e t = (pure ICM plasma) to avoid any contamination by AGN plasma. Near to and 
within the jet launching region, < lj e t, this condition was, of course, not met while the jets were 
active. For this reason there are no data points at small radii in Figures [5] and [6] when the jets were 
active. 

The initial conditions in the upper left panel of Figure [5] reveal a bias towards lower density 
within the inner core radius (< 50 kpc) with a fractional error e p ~ 15%. This is a result of the 



bias toward higher temperatures in the determination of Tjcm,i discussed in £4.1.1 By holding 
j u constant in Equation [7] for an observation with hvi oca i approximately equal to the actual gas 
temperature it can be shown that an overestimate of the gas temperature will result in an under- 
estimate of the electron density. After jet activity terminated in the RE simulation, the observed 
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density profile matches the actual ICM profile with an error, e p ~ 30%. The largest contribution 
to the error is from regions influenced by the bow shock. Excluding these regions gives e p < 15%. 
Evidence of the bow shock is seen in the actual density profile as a bump in excess of the smooth, 
observationally determined profile at r p = 35 kpc for 26.3 Myr, 70 kpc for 52.5 Myr, and 200 kpc at 
157.5 Myr. These distances correspond to the cluster-centric distance of the bow shock orthogonal 
to the jet axis. In general, the observed distribution closely matches the actual distribution for 
radii outside the inner core radius, Rci, at all times also with e p < 10%. 

The 113 profiles in Figure [6j similarly display fractional errors for the observed distribution 
e p < 15% at all times exterior to Rci- Evidence of the bow shock is seen here as well in the actual 
profile as a bump in excess over the observed profile at b = 20 kpc for 26.3 Myr, 55 kpc for 52.5 
Myr, 100 kpc for 105 Myr, and 200 kpc for 170.6 Myr. The electron density profile of the ICM 
obtained from the brightness profile was reliable to within ~ 20% outside of regions influenced by 
shocks regardless of jet intermittency and observed inclination. Inside of shock influenced regions 
the fractional error was as high as ~ 40%. 



4-1-3. Cluster Pressure Profile 

The (azimuthally averaged) radial ICM pressure profile was calculated for each observation 
from the double /3-profile model temperature and density profiles just outlined, assuming an ideal 
gas equation of state. Figures [7] and [8] show example pressure profiles determined from the observed 
temperature and density profiles along with azimuthally averaged ICM pressures and AGN pressures 
extracted directly from the RE and 113 simulation. Following the procedure used for density, only 
zones with Cj e t = (pure ICM plasma) were used to measure the ICM pressure profile. Zones 
with Cj e t > 0.01 were used to measure the AGN pressure profile. The top left panel of Figure [7] 
shows the results of the double /3-profile inversion (the "observed" profile) for the initial conditions. 
Within Rci the observed profile underestimates the actual ICM profile by > 10%. This is due to 



the underestimate of the electron density discussed in £4.1.2 Outside of Rci the ICM pressure 
profile is measured with ep < 10%. At the time the jets are turned off for the RE simulation, the 
top right panel of Figure [7| the ICM pressure is measured with ep < 20% at all radii. The signature 
of the bow shock in the ICM (solid line) can be seen at 35 kpc as a bump in excess over the observed 
profile. At 52.5 Myr, the lower left panel, the observed profile significantly overestimates the ICM 
pressure by ep > 10% within Rci- The pressure ~ 35 kpc behind the bow shock has dropped ~ 
20% from the initial conditions at this time as seen in Figure [9j By 157.5 Myr, the ICM has relaxed 
closer to equilibrium, and the observed profile measures the actual pressure to ep < 15%. 

The observed pressure profiles for the 113 simulation shown in Figure [8] reproduced the ICM 
profile to ep < 45% at all times, and at distances > 20 — 30 kpc from cluster center the observed 
profile was typically much better than that. The intermittency of the jets in the 113 run produced a 
more complex pressure distribution than the RE simulation. It cannot be captured by the smooth 
profile produced by the double f3— profile inversion. At 52.5 and 105 Myr, the actual pressure varies 
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±30% from the observed pressure within Rc2- By 170.6 Myr, the strength of the bow shock has 
diminished, and the error on the observed profile falls to ep < 10% outside Rci- Inside of Rci-, 
however, the effects of the AGN activity produces a pressure structure poorly reproduced by the 
observed profile. The measurement of the ICM pressure profile from observation was reliable to 
within ~ 20% outside of regions strongly affected by jet related shocks. Inside of shocked regions 
the measurement was only reliable to within ~ 60%. This was true for both the RE and 113 
simulations regardless of the observation orientation. 

Following convention, the observed ICM pressure profiles were used to calculate X-ray cavity 
enthalpy on the assumption that the ICM and cavity pressures were equal. The dotted lines in 
Figures [7] and [8] show the average pressure in AGN plasma at each radial bin. This pressure could 
only be observationally measured if the cavity were in exact pressure balance with the ICM and 
the exact cluster-centric distance of the cavity was known. Here we discuss how closely observation 



matches the AGN plasma profile assuming the cavity location is known. £4.1.5 discusses the effect 
of projection on inferred pressure. 

The AGN plasma pressure in the RE simulation, as shown in Figure [7j roughly follows the ICM 
pressure at 26.3 and 157.5 Myr except for high pressure at the ends of the jets where momentum 
flows drive the cavities outward (OJ10). At 52.5 Myr, the AGN plasma pressure differs by as 
much as a factor of three from the actual ICM pressure from 30-100 kpc. This discrepancy can be 
explained by the influence of the bow shock. Referring to Figure [9| shocked ICM material between 
projected distances of 30-100 kpc raised the average ICM pressure over the lower pressure inside 
of the jet cocoon. The observed pressure profile, which is sensitive to the ICM pressure, is ~ 75% 
greater than the AGN plasma pressure within Rci at this time. 

For the 113 simulation at 26.3 Myr in Figure [8j there is a significant difference between the 
AGN and ICM profiles from 25-35 kpc also due to the effects of the bow shock. At this time the 
observed pressure profile overestimates the AGN plasma pressure by ~ 50% within 20 kpc. At 52.5 
Myr, the AGN and ICM profiles approximately agree with the exception at the ends of the jets. 
Here the observed profile matches the actual AGN plasma pressure to within 13%. The observed, 
AGN, and ICM profiles all agree at 52.5 Myr to within ~ 30% except the ends of the jets. The 
intermittency of the jets impacts how well the observed profile reproduces the AGN pressure. At 
105 Myr, when the jets are inactive, the observed pressure overestimates the AGN pressure within 
30 kpc while at 170.6 Myr, when the jets are active, it underestimates it. Inferring AGN plasma 
pressure from observational measurement of the ICM pressure at a specific radius was reliable only 
to ~ 75% from the RE and 113 observations. 



4- 1-4- Cavity System Volume 



Cavity volumes were estimated from 1.5-2.5 keV observations at each analyzed epoch for both 
RE and 113 simulations. As already noted, each cavity was fit by eye with a set of ellipses (Birzan 
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et al. 2004). For each projection angle, i = 80 ,45 , and 30 , the cavities were assumed in this 
measurement to be in the plane of the sky in order to test the effects of an inaccurate or unknown 
value for the inclination. These observed cavity volumes were calculated assuming them to be 
ellipsoids of revolution around the major axis of each projected ellipse. By using multiple ellipses 
to cover each projected cavity we were better able to define the outer edge of the X-ray cavity. 



Figure 10 shows an example of the area enclosed by the ellipses chosen for the RE simulation at 
131.3 Myr observed at i = 80°. In general, it was difficult to define the edge of the cavities for the 
RE simulation once the cavities extended past Rci- For the 113 simulation the cavities were often 
outlined with a bright rim (see Figure [2]), making the edge (taken to be the inside of the rim) easier 
to find. 

The actual cavity volumes were computed by integrating the volume in the simulation data 
with Cj e t > 0.01 (partially AGN plasma). Observed and actual volumes are compared in Figure 
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Since we expect a projection bias due to foreshortening along the jet axis (see AppendixlB]), we plot 



the observed volume divided by a p /a, where a p is given by Equation B6 and a is the actual length of 
a single best fit ellipse, normalized by the actual volume. The scatter in the measurements without 
correcting for this projection bias was ~ 50%. Two features stand out in the comparisons in Figure 
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First, the a p /a correction reduced the scatter due to the projection bias to approximately 
10-15% for both RE and 113. The second obvious feature of the comparison is that the observed 
cavity volume estimates tend to be modestly smaller than the actual volumes for RE but not for 
113. The reason for this has to due with the different shapes of the cavities between RE and 
113. 113 retains a nearly elliptical area at all times while RE developed a non-elliptical shape (see 
Figures [T] and [2]) , which required many ellipses to fit. Fitting ellipses will tend to underestimate 
a non-elliptical shape if the observer requires that none of the fits extend beyond the cavity edge. 
Despite this limitation, and the subjective, observer-dependent nature of the process, our observed 
volume estimates generally agree with the actual volumes to within about ±50% (omitting the a p /a 
correction). 

The observations used to determine cavity analysis did not include any noise representing 
X-ray counts or intrumental effects. Low counts at large cluster centered distances would make 
cavity edges in those regions more difficult to identify. The long axis of all of the observed cavities, 
roughly aligned with the jet, would likely be underestimated given these conditions, which would 
reduce the measured volume by the same amount. 



4-1.5. Cavity System Enthalpy 

The observational estimates for the ICM pressure profile and cavity system volume presented 
above were used to derive the total cavity enthalpy, H b s . The cavity volume enclosed by the 
chosen ellipses for each observation was discretized into 1 kpc 3 volume elements corresponding to 
the 1 arc sec resolution of the images for a cluster distance, Dl = 240 Mpc. The pressure in each 
cavity volume was determined from the observed pressure profile and the projected cluster-centric 
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distance of the volume elements. The total enthalpy in the cavity system for each observation was 
given by 



r -^elements 

H o*s =2 S PiO-kP<?)- (16) 

i 

Figure|i~2|shows comparisons of the total energy added to the simulation volumes by jets, AEtot, 
the actual enthalpy in the cavity systems, H act , and the above observed cavity enthalpy estimates, 
H b s , throughout the RE and 113 simulations at three inclinations. The value for H ac t was computed 
from the pressure and total volume of each voxel in the simulation data with Cj e t > 0.01 to be 
consistent with the synthetic observations (see ^3]). This conservative cut meant that some ICM 
enthalpy was included in H act values, making it possible for H act > AE to t- 

Several features stand out in comparing the various energy measures. The first feature is that 
H b s , H ac t, and AEtot all agree with each other within about a factor of two for a given simulation 
and inclination angle. Second, the comparisons of enthalpy, both observed and actual, and the 
total energy measures is better for the intermittent jet simulation, 113, than for the terminated, 
RE, case. This is consistent with the analysis of the simulations reported in OJ10. In particular, 
they noted that about 50% of the jet energy (AE to t) injected during the active phase of either 
the 113 or RE model had been converted into ICM thermal or kinetic energy by the end of the 
simulation. The remaining energy increment was mostly gravitational potential energy in the ICM 
or thermal energy in the jet cocoorj^J which is roughly what H b s measures. Approximately 30% of 
AE^t ended up as gravitational potential energy in the ICM for the RE simulation. By contrast, 
a much smaller fraction, ~ 15%, of AEtot in the intermittent jet, 113, simulation is converted to 
ICM gravitational by the simulation's end. Thus, we should expect a closer match between AEtot 
and H ac t in that case. 



Another striking feature for both panels of Figure 12 is the consistency in H b s among the 
inclination angles for a given simulation and epoch. This is due to two competing projections 
effects. In particular, the estimated volume generally decreases with decreasing inclination angle 
as discussed in § |4.1.4[ while the projected distance from cluster center decreases, making a cavity 
appear to be in a higher pressure environment then it actually is. To see the net effect refer to 
Figures [7] and [8| It is evident that the ICM pressure profile can be approximated as a power 
law, P = Po(r/ro) a , over the projected distances the cavities occupy (40-200 kpc) for most of 
the simulated time. Then, assuming a given cavity is a cylinder with radius R extending from a 
projected distance r\ sin i to T2 sin i (see Appendix IB] regarding the use of sin i) and estimating 



Relatively smaller energy increments are contained at a given time in jet kinetic energy and magnetic fields. 
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a ~ -1 we would predict that the enthalpy H would be 



H = TT-R 2 p 



T2 sin i 



r\ sin i 
>'2 



-1 



dr 



tt^-R 2 p ln 
2 \r\ 



(17) 
(18) 



which is independent of i. 



For the RE simulation on the left panel of Figure 12 H Q b s always underestimated H ac t while 
the jets were active. Recall that H Q b s depends on the observed estimate of the ICM pressure 
distribution and, from §4.1, 3[ that the observed pressure profile underestimated the AGN plasma 
pressure (the cavity pressure) while the jets were active. In short, during those times the cavities 
are over- pressured as they drive moderate strength shocks into the ICM. This is consistent with 
comparisons shown in Figure [7} Further into the simulation we see H ac t declining. The cavities 
are rising buoyantly in the cluster while maintaining approximate pressure equilibrium (see OJ10). 
The thermal energy in the cavities dropped as this energy was transferred to gravitational potential 
energy. The observed values follow this trend, remaining within ~ 50% H ac t. 



The 113 simulation on the right panel of Figure 12 shows a different evolution of H ac t and H b. s 
because of the different AGN history. There is a step-like growth of H ac t due to the intermittency of 
the jets. The time delay between the peaks of H ac t and H b s at each step is due to the difference in 
evolution of the observed ICM pressure and the actual AGN plasma profiles. When the jets turn on 
the cavities are over-pressured with respect to the ICM, but they eventually expand to approximate 
pressure balance. Prior to the expansion of the cavities, however, the higher energy content of the 



cavities cannot be accurately measured by the procedure described in £4.1.3 Therefore, an increase 
of H bs lags behind an increase of H act . 



4.2. Ages 



Three characteristic timescales are commonly employed for determining cavity age. For a 
cavity centered at a projected distance r p from cluster center, radius R, cross section S, drag 
coefficient C, and volume V these times are: 1) the buoyant rise time ty 



■uoy 



1 ^CS/{2gV), 2 



the "refill time" t r = 2y/R/g, and 3) the sound crossing time t c = r p /c s ( Birzan et qZ.||2004 ) . These 
times can be compared to known ages from the simulations given measurements of the sound speed, 
c s , and g from the synthetic observations. Measurements of each age were made from observations 
of the N and S cavities (see Figures [I] and [2J at the end of the RE simulation, t = 157.5 Myr, 
and 113 simulation, t = 170.6 Myr, at three different inclination angles. In these observations we 
represented each cavity as a single ellipsoid with semi-major axis a and semi-minor axis b. Following 



a procedure similar to Birzan et al. (2004), the value of r p was the projected distance from cluster 
center to the center of the cavity, the radius was given by R = Vab, and the cross section was given 



by S 



^max' where 



was half the maximum azimuthal width of the cavity. The volume V 
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was determined for each cavity following the method described in £4.1.4 The sound speed was 
given by c s = \p~jk < Tjcm > UJ^h)- At the end of the RE simulation < Tjcm >= 2.77 keV 
giving c s = 941 km s , and at the end of the 113 simulation < Tjcm >= 2.82 keV, giving c s = 
950 km s _1 . Refer to Appendix [c] for details on estimating g, which we assumed to be constant, 
from the synthetic observations. We let C = 1 for simplicity. 

Table [l] shows the measured parameters and age estimates for each observation at the end 
of the RE and 113 simulations. Both tbuoy and t c are affected by projection. For this reason, 
we would expect the buoyant rise time to vary as £& uo2/ oc r p /y/a oc y/sini due to projection 
effects on both the projected distance and semi-major axis of the cavity. The sound crossing 
time, however, should vary more rapidly with inclination as t c oc r p oc sini. The cavity pairs 
for both simulations approximately show this trend for tbuoy and t c . The refill time shows a 
weaker dependence on inclination angle because t r oc a/ (aft) 1 / 2 oc (sini) 1 ^ (recall that we assumed 
constant g in calculating ages). An important aspect of Table [l] was that for most cases tbuoy < t c , 
which implies a terminal buoyant velocity greater than the sound speed. This unphysical result 
could have been avoided in a number of ways. In their analysis of Hydra A, for example, Wise et 



al. (2007) represented the cavity system as a series of spherical bubbles. Approximating the end 



of the N cavity of the RE simulation observed at i = 80° as a sphere would increase r p to ~ 200 
kpc and would decrease V to (4/3)7rr 3 ~ 1.4 x 10 5 kpc 3 . Given these measurements for an outer 
spherical cavity, tbuoy ~ 244 Myr while t c ~ 208 Myr. Another approach may have been to assume 
C > 1 similar to values empirically estimated by Jones & DeYoung (2005). The parameters and 
measurements used in a buoyant rise model are not very well constrained. A buoyant model also 
did not properly capture the evolution of the cavities in the RE or 113 simulations (see OJ10). For 
these reasons we chose not to use tbuoy as the cavity age. The equations for t r and t c are related, 
and the models only differ in the length over which material moves. By convention, t r uses the 
cavity radius, which is not affected by projection. We instead use the values of t c for the cavity 
ages in subsequent calculations so that our analysis demonstrates the dependence on projection. 
When projection did not greatly effect our measurements t c was reliable to within ±20%. 



4.3. Cavity Power 



The enthalpy and age of a cavity system are typically combined into a characteristic power 
called the cavity power P C av — H/t. If Pcav can be deposited into the ICM it may balance the 
cooling in the host cluster. It is therefore common practice to compare Pcav with, the cooling [e.g., 
Birzan et qfl2004| |Rafferty et of.|2006| |Cavagnolo et aZ.|[2010| ). The cavity power should represent 
a lower limit to the total luminosity of the AGN because it does not account for energy already 
deposited into the ICM or other forms of AGN plasma energy such as magnetic, potential, or kinetic 
(see j |4.1.5D . Table g shows a comparison of the observed P ca v,obs with the actual total average jet 
luminosity, Lj e t t act, for both RE and 113 observed at the three different inclinations at the end of 
each simulation. P cav ,obs was computed as H b s /(t), where (t) was the average of the ages in Table 
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[Tjfor the N and S cavities, while the mean jet power, Lj e t :ac t = AE(t)/t, where t = 157.5 Myr and 
t = 170.6 Myr for the RE and 113 simulations respectively. The differences between P ca v,obs and 
Lj e t,sim at these times were characteristic of early times. 

Observations of the RE and 113 simulations produced P ca v,obs increasing with decreasing in- 



clination primarily due to projection effects on measuring the ages as discussed in {4.2 For the 
RE simulation, P cav ,obs < Lj et ^ct as we would expect for all but the smallest inclination angles. 
The 113 observations, however, resulted in P cav ,obs > Ljet,act at all orientations. The underestimate 
of cavity system age was the dominant reason P ca v,obs exceeded Lj e t t act in these observations. For 
measurements from observations at i = 80°, when the error on the assumed inclination is small, 
Pcav,obs w as within ~ 40% of Lj etA ct, and P ca v,obs was within a factor of three of Lj etA ct across all 
of the observed inclination angles. 



5. Conclusions 

We have presented an analysis of the reliability of common techniques used to extract X-ray 
cavity enthalpy, age, and mechanical luminosity from X-ray observations of cavity systems. By 
utilizing synthetic X-ray observations of detailed simulations we were able to directly compare 
observationally determined and actual values from the simulations. The important results from 
this work are: 

• The synthetic observations of the 113 simulation show bright rims outlining the cavities at 
each analyzed epoch out to 170.6 Myr while the RE simulation does not show bright rims. The 
difference in the AGN history represented by each model accounts for this difference. 113 had 
periodic injection of energy into the cavities throughout the simulation while RE deposited all of 
its energy early on. 

• Observationally measuring X-ray cavity enthalpy is reliable to within approximately a factor 
of two across a wide range of age and inclinations for the models of jet intermittency presented 
here. Several steps go into determining the enthalpy, and each may introduce significant errors. 
Extracting the ICM electron density profile, for example, was reliable to within ~ 20% outside 
of regions strongly influenced by shocks. Inside recently shock influenced regions the error was ~ 
40%. Combining this and the temperature profile into the pressure inside of the cavity at a given 
cluster-centric distance may not be as accurate. During periods of jet activity, the observationally 
determined pressure may differ by as much as ~ 75% from the cavity pressure. This is related 
to the supersonic speeds of the jets through the ICM and the consequent post-shock pressure 
enhancement. Our measurements of cavity volume were within ±50% of the actual total cavity 
system volume. This process is subjective, however, and a more robust and objective method for 
finding and outlining cavities should be developed. The overall effect of each of these measurements 
is contained in the factor of two reliability of enthalpy measurements. The energy required to offset 
cooling in clusters can be characterized as r/PV. An approximate factor of two span in r/ in our 
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tests is due largely to uncertainties in the measurement of PV. 

• The determination of cavity age from one or more of the commonly used age estimates could 
potentially be misleading. The buoyant rise model was not an accurate description of the evolution 
of the cavities in our simulations. A simple application of this model implied unrealistic terminal 
velocities greater than the sound speed. The refill time model produced ages within ±15% of the 
correct age regardless of the error on assumed inclination. It relied on an accurate measurement of 
the gravitational acceleration, however, which assumes the cluster to be in hydrostatic equilibrium. 
This assumption may be not be valid for a given cluster. We preferred to use the sound crossing 
time as a simple and fairly robust model for cavity age. For a well constrained inclination angle, 
our measurements were within ±20% of the actual cavity age. 

• Observationally measuring the cavity power produced values within a factor of ~ 3 of the 
average total jet luminosity from our simulations regardless of assumed inclination angle. The 
observed cavity power was within 40% of the jet luminosity if the projection effects were negligible. 
At all observed inclination angles for the 113 simulation and i = 30° for the RE simulation the 
cavity power overestimates the average jet luminosity largely due to underestimates in the cavity 
system age. 

This work was supported at the University of Minnesota by NSF grant AST0908668 and by 
the University of Minnesota Supercomputing Institute. PJM and SMO were supported in part by 
the Graduate Dissertation Fellowship at the University of Minnesota. SMO was also supported by 
NASA Astrophysics Theory Program Grant NNX09AG02G. We are grateful to Brian McNamara 
and Paul Nulsen for very fruitful conversations and to an anonymous referee for help in improving 
the original manuscript. 



A. Fitting (3 Profiles 

To calculate a best fit to the brightness distributions comparable with observations, the flux 
from the synthetic observation was scaled by the net counts from Chandra observations of Hydra. 
Counts were taken from the evt2 files in the Chandra archive (ObsIDs 575 and 576, Chandra 
calibration program; ObsIDs 4969 and 4970, program 05800556, P.I. McNamara) ( |Wise et al. 



2007). The total number of counts with energies from 1.5-2.5 keV in a circular region centered 
on the cluster with a diameter equal to the size of the synthetic observations, 12.3 arc min a side, 
was extracted for each observation. A total of N c = 451345 counts for the four observations were 
detected in this region over 227 ksec. The scaling factor £ was calculated as 

N c counts 

1 * pixels 

Si erg s -1 cm~ 2 
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where Si is the flux for a given pixel observed over a band from 1.5 to 2.5 keV. Xnorm was then 
calculated from the simulated average counts in a bin, the product of £{Si), assuming Poisson 
statistics. The downhill simplex method was used to minimize Xnorm by simultaneously fitting 
Soi, Rc,i, and Both the RE and 113 grids extended only 300 kpc. To ensure that our fits 
were not affected by the truncation of the atmosphere, additional azimuthally averaged data points 
were generated from a synthetic observation of the cluster discussed in §2,2| extending the initial 
conditions to 600 kpc. 



B. Projected Semi-major Axis of an Ellipsoid 



The projected semi-major axis of revolution of an ellipsoid can be determined from the incli- 
nation angle, i, and the un-projected aspect ratio taken to be ip = a/b, where a is the semi-major 
axis and b is the semi-minor axis along the other two dimensions. Figure |i~3| shows an example of 
the geometry of the problem where the line segment OC is the projected semi-major axis. We start 
with an ellipse defined by 



x 2 y 2 
a 2+ ¥ 



1 



(Bl) 



lying in a plane that includes the line of sight with a tangent line of slope 
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from the observer. The coordinates D = {xo,yo}, where the tangent line intersects the ellipse, are 
given by 
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The equation for the tangent line is then 

y = y + tani(x - x ). 

This line will intersect the orthogonal line passing through the origin given by 

x 

y 



tan i 



Solving for x and y where Equations B4 and B5 intersect one finds the length OC to be 



OC 



h + tan 2 i) (l + tan 2 i 
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Figure 14 shows the projected semi-major axis normalized by the actual semi-major axis, a p /a, 



as a function of i for an ip = 2 and tp = 4. At large inclination angles (i > 45°) sini approximates 
a p /a to within 12%. For simplicity we typically chose to correct an observed semi-major axis in 
that regime with sini, as it does not require that ip be known. 
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C. Gravitational Acceleration 



The gravitational acceleration at a given projected cluster radius r p can be determined from 
the synthetic observations using an "observed" measurement of the distribution of gravitating mass. 
An estimate of the mass within a radius r p from the synthetic observations was be made assuming 
the ICM plasma was in hydrostatic equilibrium with an electron density profile derived from a 



double /3-profile (Equation 15), and dominated by the pressure of the thermal gas. ^From these 



assumptions the distribution, given by Xue & Wu (2000), is 



M(r p ) = Y,Mi 



n e (0) 



1 ■ ( 1 ^ r P+ R C, 
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The gravitational acceleration, g, was then computed as g(r p ) = GM(r p ) /r p . The resulting value 
of g varied by approximately a factor of two over most of the observed cluster (10-300 kpc). In 



deriving the timescales discussed in ^4.2 a single value for g is required. We adopted the value of g 



to be a simple average over the observed g{r p ) profile. This value for g was within a factor two of 
the actual gravitational acceleration used in the simulation, g a ct( r p) (Equation [21), over the range 
10 kpc <r p < 300 kpc. 
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Fig. 1. — Synthetic observations of the RE simulation at three epochs and three orientations 
from the energy band 1.5-2.5 keV. Each row corresponds to one inclination angle with the time 
progressing from left to right. The observations were divided by a best-fit double (3— profile to 
emphasize the X-ray cavities. The N and S cavity labels distinguish the cavities from each jet. 
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Fig. 2. — Synthetic observations as in Figure [T] of the 113 simulation. 
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Fig. 3. — Double f3— profile fits to azimuthally averaged radial brightness brightness distributions 
measured from the synthetic X-ray observations. On the left is the fit for the initialized cluster 
of both RE and 113 observed at an inclination of 80°. On the right is a fit for the RE simulation 
at 157 Myr observed at an inclination of 45°. The bottom panel of each plot shows the fractional 
error of the fit to the observed distribution. 
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Fig. 4. — The observed and actual temperature profiles from the initialized cluster of both the RE 
and 113 simulations. The bottom shows the fractional error of the observationally derived profile 
from the actual profile. 
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Fig. 5. — ICM electron density profile for the RE simulation at four characteristic epochs at an 
inclination of 45°. The solid line shows the actual azimuthal average measured from the simulation. 
The dash- dot line is the profile determined from the inverted double /3-profile fits to the synthetic 
X-ray observations. The bottom of each panel shows the fractional error of the observed profile 
from the actual profile. 




Fig. 6.— 



Similar to Figure [5] for the 113 simulation. 
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Fig. 7. — ICM pressure profile for the RE simulation at four characteristic epochs at an inclination 
of 45°. The solid line shows the actual azimuthal average measured from the simulation. The 
dash-dot line is the profile determined from the inverted double /3-profile fits to the synthetic X-ray 
observations. The dotted line is the average pressure in AGN plasma at that radius measured from 
the simulation. The bottom of each panel shows the fractional error of the observed profile to the 
actual ICM profile. 




Fig. 8.— 



Similar to Figure [7] for the 113 simulation. 
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52.5 Myr 




Fig. 9. — A slice through the midplane of the computational box showing the ratio of the pressure 
at 52.5 Myr into the RE simulation to the initial conditions. The pressure within 50 kpc of the 
cluster center has dropped by ~20% from the initial conditions. The bow shock is visible as an 
increase in pressure by a factor of ~ 2.5 from the initial conditions. (A color version of this figure 
is available in the online journal.) 
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Fig. 10. — Example of area enclosed by ellipses fit by eye (solid line) to a 1.5-2.5 keV observation 
of the RE simulation at 131.3 Myr divided by a best-fit double (3— profile. (A color version of this 
figure is available in the online journal.) 
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Fig. 11. — The observationally determined total cavity volume to the actual total cavity volumes 
from each every analyzed epoch for the RE simulation (left) and the 113 simulation (right). (A 
color version of this figure is available in the online journal.) 
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Fig. 12. — Comparison between the actual enthalpy in AGN plasma, i^actj an d observed enthalpy, 
-ffobsj for the RE simulation (left) and the 113 simulation (right). H D b s measured from observations 
at three inclination angles are shown. The total energy added to the computational grid by the 
AGN, AEtot, is shown as a dotted line for reference. 



-33- 



Table 1. Cavity Ages 



Model 


i 


Cavity 


r P 


a 


b 


bmax 


V 


9 


^buoy 




t c 








(kpc) 


(kpc) 


(kpc) 


(kpc) 


(10 5 kpc 3 ) 


(10~ 8 cm s~ 2 ) 


(Myr) 


(Myr) 


(Myr) 


RE 


80° 


N 


146 


118 


25 


32 


3.69 


2.4 


110 


167 


152 


RE 


80° 


S 


134 


115 


26 


34 


3.62 


2.4 


108 


168 


139 


RE 


45° 


N 


97 


89 


27 


32 


2.43 


2.6 


86 


153 


101 


RE 


45° 


S 


97 


87 


27 


35 


2.58 


2.6 


91 


152 


101 


RE 


30° 


N 


74 


70 


24 


33 


2.09 


2.6 


73 


140 


77 


RE 


30° 


S 


74 


68 


24 


38 


2.06 


2.6 


85 


139 


77 


113 


80° 


N 


135 


120 


42 


43 


9.95 


2.4 


83 


191 


139 


CO 


80° 


S 


138 


120 


45 


46 


9.83 


2.4 


91 


195 


142 


113 


45° 


N 


96 


88 


40 


43 


9.02 


2.5 


61 


172 


99 


113 


45° 


S 


95 


90 


43 


43 


8.56 


2.5 


62 


176 


99 


113 


30° 


N 


69 


69 


41 


42 


6.63 


2.5 


50 


162 


72 


113 


30° 


S 


73 


71 


41 


42 


6.67 


2.5 


52 


164 


76 



Note. — Several cavity age estimates measured from observation for both RE and 113 at 157.5 Myr and 170.6 Myr 
respectively. 



Table 2. Cavity Power 



Model 


i 


H 

(10 60 erg) 


(t) 
(Myr) 


-P 'cav,obs 

(10 44 erg s- 1 ) 


(10 44 erg s- 1 ) 


RE 


80° 


3.3 


146 


7.2 


12 


RE 


45° 


3.4 


101 


10.6 


12 


RE 


30° 


3.6 


77 


14.9 


12 


113 


80° 


11.6 


141 


26.1 


21 


113 


45° 


14.4 


99 


46.1 


21 


113 


30° 


13.9 


74 


59.5 


21 



Note. — Cavity power P cav ,obs measured from observation for both 
RE and 113 at 157.5 Myr and 170.6 Myr respectively. The average 
jet luminosity in the simulation out to those respective ages, Lj et;act , 
is given for comparison. 
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Fig. 13. — Diagram showing an ellipse with an aspect ratio of 4 observed at an inclination angle 

O - 

i = 45 . The length of the line segment from point OC is the projected semi-major axis. The 
segment is defined by a tangent line and the orthogonal line passing through the origin. 
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Fig. 14. — The projected semi-major axis of an ellipse with an aspect ratio of 2 (solid curve) and 
aspect ratio of 4 (dashed curve) as a function of inclination angle. A sini (dash-dot curve) curve 
is shown for comparison. At large inclination angles a simple sin curve closely approximates the 
projected semi-major axis for both aspect ratios. 



